Astronomy & Astrophysics manuscript no. postagb 


© ESO 2008 


February 5, 2008 





A 1.2 mm MAMBO survey of Post-AGB stars 

C. S. Buemi 1 , G. Umana 1 , C. Trigilio 1 , and P. Leto 2 

INAF, Osservatorio Astrofisico di Catania via S. Sofia, 78 95123, Catania Italy 
2 INAF, Istituto di Radioastronomia Sezione di Noto, CP. 161 96017, Noto (SR) Italy 



Received 

ABSTRACT 

Aims. We performed a millimetric survey of a sample of 24 post-AGB stars aimed to search for emission from circumstellar 
matter, in order to investigate the physical properties of the outer parts of the envelopes. 

Methods. The observations were conducted using the 37-channel Max-Planck Millimeter Bolometer array at the 30-meter 
IRAM telescope. The continuum emission toward the detected sources was used to quantify the mass of the emitting dust. We 
combined our observations with data available in literature to construct the spectral energy distribution (SED) of the sources. 
When the observational data cover a spectral range wide enough, some properties of circumstellar envelopes were derived by 
comparison with spectra computed using a radiative transfer code. 

Results. Of the 24 objects in our sample, we detected millimetric continuum emission toward 11 sources. Two other sources were 
detected at flux level close to 3a. The derived circumstellar dust masses range between 0.4 and 24 x 10 -4 Mq, but these results 
are affected from the uncertain about the source distances. The parameters derived from the SED fits are consistent with the 
values characteristic for these kind of object. As confirmed from the flux density extrapolated in the first light channels of the 
Atacama Large Millimetric Array, such sources would be good targets for future high resolution mapping with the ALMA facility. 
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1. Introduction 

Stars on the red giant branch undergo a process of copious 
mass-loss leading to the formation of a circumstellar en- 
velope (CSE). These envelopes, quilted up by the ejected 
material, contain atomic and molecular gas and are char- 
acteristic of the post Asymptotic Giant Branch (AGB) 
evolutionary phase, which results in the formation of a 
Planetary Nebula (PNe) . Because of the huge amount of 
processed material returned to the ISM, this evolutionary 
phase is very important for the chemical evolution of the 
Galaxy. Yet, the short transition phase between the end 
of the AGB and the formation of a new PNe is still poorly 
understood. In particular, it is quite challenging to under- 
stand how the almost symmetric CSE observed around 
AGB stars transform themselves in the highly structured 
morphologies observed in high dynamical range optical 
images of PNs. 

The importance of studying objects in the way to PNe 
resides in the fact that in their CSE the unknown physical 
mechanisms, that shape the PNe, are already at work, as 
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HST images of multi-polar outflows in post-AGB stars 
appear to indicate (Sahai[2001 ). 

Despite of numerous efforts to identify the shaping 
agent/s in PNs (Sahai & Trauger 19981 Garcia-Segura et 
al. l2005fl . there are no observative evidences that support 
one mechanism in favors to the others (Balick and Franck 
(TUU2|) . 

Dust is ubiquitous detected in post-AGB and PNe and, 
quite often, it is in the form of disks and tori. Some models 
evocate these structures, whose origin is however under 
debate, as important ingredient for the mechanism which 
produces collimated outflows observed in the CSE (Frank 
et al. 1997; Huggins et al. I2004J1 . In this contest, very 
important informations can be provided by studying the 
physical properties of CSE and in particular, the presence 
and the spatial distribution of the circumstellar dust, as 
this may help in finding clues to understand the nature of 
the shaping agent. 

Detailed studies of the spectral energy distribution 
(SED) of a small sample of optically bright post-AGB can- 
didate stars have shown that these objects can be divided 
into two groups, depending on the shape of the IR excess 
(Trams et al. 119911 van der Veen et al. H994fl : sources with 
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Table 1. Properties of stars in our sample. IRAS fluxes come directly from Point Source Catalog (PSC). 



Sources with Optical Counterpart 



IRAS Name 


RA(J2000) 


DEC(J2000) 


Fia 


F 25 


F 6 o 


Fioo 








[Jy] 


[Jy] 


[Jy] 


[Jy] 


00210+6221 


00 23 51.2 


+62:38:16 


48.5 


51.9 


12.5 


<23.2 


01174+6110 


01:20:44.9 


+61:26:18 


4.1 


16.9 


33.9 


4.1 


04296+3429 


04:32:56.6 


+34:36:11 


12.7 


45.9 


15.4 


<9.2 


05089+0459 


05:11:36.1 


+05:03:26 


7.4 


21.9 


11.9 


3.8 


06530-0213 


06:55:32.1 


-02:17:30 


6.1 


27.4 


15.1 


4.1 


07134+1005 


07:16:10.2 


+09:59:48 


24.5 


116.7 


50.1 


18.7 


07331+0021 


07:35:41.1 


+00:14:58 


15.3 


68.1 


18.5 


3.7 


17436+5003 


17:44:55.4 


+50:02:39 


6.1 


184.0 


152.0 


48.7 


19114+0002 


19:13:58.6 


+00:07:31 


31.3 


648.3 


515.9 


168.1 


20000+3239 


20:01:59.4 


+32:47:32 


15.0 


71.0 


30.0 


<43.1 


20028+3910 


20:04:35.0 


+39:18:38 


41.8 


210.8 


143.1 


46.5 


22223+4327 


22:24:30.6 


+43:43:03 


2.1 


37.1 


22.4 


9.5 


22272+5435 


22:29:10.3 


+54:51:06 


73.9 


302.4 


96.6 


41.0 


22574+6609 


22:59:18.4 


+66:25:48 


9.0 


29.5 


10.7 


2.5 


23304+6147 


23:32:45.0 


+62:03:49 


11.4 


59.1 


26.6 


7.2 



Sources without Optical Counterpart 



07430+1115 


07:45:49.8 


+ 11:08:25 


7.7 


22.9 


10.7 


2.5 


18454+0001 


18:48:01.5 


+00:04:47 


10.8 


14.5 


13.6 


<384 


18514+0019 


18:53:57.9 


+00:23:24 


4.9 


23.4 


17.3 


<152 


18576+0341 


19:00:11.2 


+03:45:46 


58.5 


425.0 


274.7 


<1660 


19024+0044 


19:05:01.5 


+00:48:48 


2.9 


48.8 


42.5 


15.7 


19075+0432 


19:10:00.0 


+04:37:06 


5.2 


28.1 


31.8 


14.4 


19454+2920 


19:47:24.3 


+29:28:12 


17.3 


89.6 


54.4 


14.7 


20144+4656 


20:15:58.3 


+47:05:39 


1.2 


17.0 


20.0 


<85 


21537+6435 


21:55:04.6 


+64:49:54 


6.9 


26.1 


13.3 


<6.1 



a broad IR excess extending from the near infrared until 
the far-IR have both hot and cool dust in their circum- 
stellar shells and sources with only a far-IR excess show 
only the presence of cool dust. The double peaks in the 
SEDs appear to be characteristic for objects in transitions, 
but the presence of only cool dust seems to point out ob- 
jects more evolved towards PNs. Thermal emission from 
dusty envelopes may extend up to millimiter wavelengths, 
depending on the temperature of the dust. In particular, 
millimetric observations are essential to better determine 
the SEDs of the sources and thus to put more stringent 
constraints to the model of their circumstellar envelopes. 
Moreover, millimetric observations will allow to assest the 
presence of multiple dusty shells to be related to different 
mass-loss episodes undergone by the stars. 

In this paper we present the results of a 1.2 mm sur- 
vey aimed to detect the millimetric emission from a sam- 
ple of post-AGB stars for a robust modeling of their cir- 
cumstellar envelopes. This appears as a first, necessary 
step in preparation for the use of the new interferometers 
(i.e. ALMA and CARMA) that, in the very near future, 
will be able to resolve, in unprecedented spatial resolution, 
the circumstellar geometry for a full understanding of the 
shaping mechanisms. 



2. Observations 

2.1. Sample selection 

In the attempt to determine a possible millimetric emis- 
sion due to thermal dust emission from the circumstellar 
envelope, we observed a sample of Post-AGB stars with 
the 30m IRAM telescope. 

Since the dusty envelopes show a clear signature in 
the far-IR spectrum of the stars, the IRAS colour-colour 
diagram has been successfully used by several authors in 
systematic searches for post-AGB objects by looking for 
sources in-between the locus of the planetary nebulae and 
late-type AGB stars. Garcia-Lario et al. I|1997[l compiled 
a list of stellar objects, characterized by strong F-IR ex- 
cess, which occupy the same region of the IRAS color dia- 
gram as AGB stars and PNs. The sample contains 126 
Post-AGB. Among those, a large fraction (70%) represents 
post-AGB in a very early stage, still heavily obscured in 
the optical. 

Another method of looking for transition objects is 
concentrated on optically bright objects with an IR ex- 
cess due to circumstellar dust (Pottash and Parthasarathy 
HUggl Trams et al. HMTl Oudmaijer et al. HM^jl . These re- 
sulted in the detection of objects scattered in the IRAS 
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Table 2. Measured millimeter fluxes. 



Sources with Optical Counterpart 



1 1 vAo in am e 


Date 


i line 


r lux 


t 

Weather 










Tel 

[s] 


[mjy] 




UUzlU+Ozzl 


21 


Jan 


02 


IzUU 


. . . ±z.4 


13 


fit 1 TA r CI in 

Ull f4+ollU 


21 


Jan 


02 


IzUU 


zz. y±i. ( 


13 


u^zyo-ro^zy 


14 


Feb 


02 


1 900 
IZUU 




A 
i i 


05089+0459 


05 


Feb 


02 


1440 


...±1.0 


A 


06530-0213 


14 


Feb 


02 


1200 


4.8+1.5 


A 


07134+1005 


07 


Feb 


02 


1200 


14.0+1.5 


A 


07331+0021 


07 


Feb 


02 


1200 


. . . ±1.9 


B 


17436+5003 


05 


Feb 


02 


1080 


15.2+1.1 


A 


19114+0002 


04 


Feb 


02 


480 


68.2+3.4 


B 


20000+3239 


21 


Jan 


02 


1200 


11.4+1.7 


B 


20028+3910 


21 


Jan 


02 


1200 


11.9+1.5 


A 


22223+4327 


04 


Feb 


02 


1200 


...±1.5 


A 


22272+5435 


04 


Feb 


02 


1200 


35.3+1.7 


A 


20144+4656 


21 


Jan 


02 


1200 


...±1.8 


B 


23304+6147 


21 


Jan 


02 


1200 


...±1.8 


B 


Sources without Optical Counterpart 


07430+1115 


07 


Feb 


02 


1200 


...±1.3 


A 


18454+0001 


22 


Jan 


02 


1200 


. . . ±1.4 


A 


18514+0019 


21 


Jan 


02 


1200 


8.2+1.3 


A 


18576+0341 


03 


Feb 


02 


840 


>45.4±2.2 


B 


19024+0044 


22 


Jan 


02 


1200 


...±1.3 


A 


19075+0432 


04 


Feb 


02 


1200 


6.2±1.2 


A 


19454+2920 


03 


Feb 


02 


1200 


. . . ±4.3 


C 


20144+4656 


05 


Feb 


02 


1200 


. . . ±2.7 


C 


21537+6435 


04 


Feb 


02 


1200 


6.3±1.6 


A 



Wheatear condition: A) good, B) poor, C) bad. 



colour-colour diagram, as they show different amounts of 
IR excesses and different IR-colours. 

Oudmaijer et al (1992 1996), looked for transition ob- 
jects with an optical counterpart by performing a cross- 
correlation of the SAO optical catalogue with the IRAS 
point-source catalogue, selecting supergiants with spectral 
type between B and G and with an IR excess due to cir- 
cumstellar dust. 

To select a sample of possible targets we used the most 
complete compilations of objects with no optical counter- 
part, from Garcia-Lario et al. (1997) and of objects with 
optical counterpart, from Oudmaijer et al (1992; 1996), 
selecting the stars identified as post-AGB, which are asso- 
ciated with highly evolved post-AGB stars with low-mass 
progenitors. We further include in the sub-sample of stars 
with optical counterpart 11 stars, originally not classified 
by Oudmjier et al. i| 1 992 ) as Post-AGB, for which there are 
compelling observational evidences that are in the post- 
AGB evolutionary stage (van der Veen et al. IT993| van 
Whicker ITM7|l . 

A small sample of post-AGB stars has been observed 
in the millimetric band with the JCMT (van der Veen 
et al. I1994|l and the detection rate appeared to be well 
correlated with the F 6 o IRAS flux, in the sense that all 



the stars with F 60 > 10 Jy have been detected. If the 
extra infrared excess has the shape of a cold black body, 
as expected from a more distant dust shell, a F 60 ~ 10 
Jy would imply a 1.3 mm flux higher than 36 mjy and 
thus easily detectable with the new Bolometer. In order 
to maximize the probability to have a detectable flux at 

1.2 mm we thus selected, from our original sample, only 
those sources with Fqo > 10 Jy. This reduces our sample 
to 34 targets. 

In this paper we report on the results relative to a 
subsamples of 24 objects, which are given in Tabled 

2.2. The 30m IRAM observations and results 

The 37-channel Max-Planck Millimeter Bolometer 
("MAMBO"; Kreysa et al. ITMSjl array at the 30-meter 
IRAM telescope on Pico Veleta (Spain) was used to 
perform the survey. The observations were made between 
21 January and 7 February 2002, using the standard 
ON-OFF technique, chopping the secondary mirror of 
the telescope by about 50" in azimuth, at a rate of 1 Hz. 
The FWHP of our beam was 10". 5 at 1.2 mm. 

For each source the observations were typically ob- 
tained in blocks of 4 or 5 scans lasting 4 minutes each. 
Frequent skydip observations have been used to determine 
atmospheric extinction as function of elevation and time. 
The data were analyzed using the MOPSI software (Zylka 
1998). The flux calibration was performed by observing 
either Mars and Uranus to determine the flux conversion 
factor. For each channel the sky noise was subtracted by 
computing the weight mean of the signals from the sur- 
rounding six channels. 

In table [21 the resulting 1.2 mm flux densities mea- 
sured for the detected objects are reported. Out of the 
24 observed sources, we carried out 11 detections over 
the threshold of 4cr and two uncertain detections at 3<r. 
Four of the detected and one of the undetected sources 
overlap the sample observed from Walmsley et al. 0991), 
by using the 30m IRAM telescope at the wavelength of 

1.3 mm, but with a different bolometer. The measured 
flux densities are in good accordance with the exception 
of IRAS 22272+5435, which shows a flux level significa- 
tively higher. 

The sample includes IRAS 01174+6110 and 
IRAS 18576+0341, which were firstly classified as 
post-AGB stars on the basis of their IRAS colours, but 
probably having different nature. IRAS 01174+6110, in 
fact, is likely an HII region (Kelly et al. 2005), while 
IRAS 18576+0341 has been recently recognized as new 
LBV (Pasquali & Comeron l2T)f]2l Clark et al. l2T)t)3l 
Umana et al. 2005). For these two objects we just report 
the measured flux density. 

In particular, IRAS 18576+0341 will be object of a 
future more detailed analysis assembled with radio high 
resolution results. It is to note that Umana et al. (2005) 
pointed out that the source's right ascension is about 4" 
shifted with respect to values previously reported in lite- 
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Fig. 1. Spectral energy distributions of the emission towards the detected sources. The asterisks indicate measures 
from this work. The dashed lines show the modified blackbody fits to the long-wavelength data, as discussed in the 
text (par. 3.f). 



rature (Garcia-Lario I1997fl . Since the first contour level 
in the 22 GHz map defines a source size close to 10", it 
is reasonable to consider the measured 1.2 mm flux as a 
lower limit. 

Some of the detected object in our sample have been 
detected in 13 CO and/or 12 CO transition J=2-l at A= 1.3 
mm, that is within the band of our observations (Hrivnak 
et al. 120051 Bujarrabal et al. 119921 Bujarrabal et al. 
I2001J) . We verified that the contribution of such emission, 
spreaded in our bandwidth of about 80 GHz, is lower than 
the errors associated to the measures for all the detected 
sources but IRAS 19114+0002 and IRAS 22272+5435. In 
particular, on the basis of the observations by Bujarrabal 
et al. ljTM3|) . we derived for IRAS 19114+0002 an emis- 
sion, in the MAMBO band, of 15 mJy from 12 CO and 
of 1.6 mJy from 13 CO; for IRAS 22272+5435 the contri- 
bution of 12 CO emission is about 7 mJy (Hrivnak et al. 
2005). In the following calculations, we thus subtracted 
these contributions from the observed fluxes. 

2.3. Spectral energy distributions 

In order to build up the spectral energy distribution 
(SED) of the detected sources, updated with our milli- 



metric data, and to investigate its corrected low frequen- 
cies shape, we combined our MAMBO observations with 
infrared (2MASS + IRAS + MSX) and optical data avail- 
able in literature. The data are corrected for the interstel- 
lar extinction, with the exception of IRAS 19075+0432 
due to the lack of informations about the extinction to- 
ward the source. Moreover we have excluded from this 
analysis IRAS 01174+6110 because of the uncertainty on 
its nature. 

For those objects for which the observations cover a 
spectral range wide enough, the resulting SEDs (fig. ^) 
show the typical double peaked shape, with the two peaks 
corresponding respectively to the optical emission from 
the photosphere of the central star and to the thermal 
emission from the circumstellar dust. 

3. Analysis 

3.1. Dust Masses 

From a theoretical point of view, the mass of the circum- 
stellar dust surrounding the central object can be derived 
from the observed continuum millimetric flux. Assuming 
that the flux density measured at this frequencies is to be 
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Table 3. Derived values for the emissivity index, dust mass, absorption coefficient and envelope dust mass. 



Source with Optical Counterpart 



IRAS Name 


T d 


P 


Xl.3 




M d 


d 




[K] 




hV 1 ] 




[M ] 


[kpc] 


07134+1005 


135 


1.11 


1.75 




8.5 (-4) 


2.4 


17436+5003 


100 


1.55 


0.88 




6.2 (-4) 


1.2 


19114+0002 


100 


1.41 


1.09 




4.3 (-2) 


6.0 


20000+3239 


140 


0.94 - 1.55 


0.88 - 2.1 


0.9 - 


2.3 (-4) 


1.0+ 


20028+3910 


100 


1.65 


0.75 




2.4 (-3) 


2.9 


22272+5435 


145 


1.03 


1.98 




6.2 (-4) 


1.6 


Source without Optical Counterpart 


18514+0019 


130 


0.89 - 2.18 


0.32 - 2.5 


0.6 - 


4.8 (-4) 


1.0+ 


19075+0432 


85 


1.37 


1.2 




1.5 (-4) 


1.0+ 


21537+6435 


140 


0.87 - 0.97 


2.18 - 2.55 


4.3 - 


5.1 (-5) 


1.0+ 



Assumed distance. 



ascribed to the thermal emission from optically thin and 
isothermal dust, the dust mass (Md) and the flux density 
are directly proportional (Hildebrand 1983 ): 

where B v (T d ) is the Planck function for dust temperature 
Td, d is the distance to the source and Xv is the dust 
opacity at the observing frequency We can thus estimate 
the dust mass by using the Rayleigh- Jeans approximation: 

M d = (1) 

The value of Xv is the major uncertainty that affects the 
conversion of the millimetric flux density in dust mass. 
Following Hildebrand's approach, we can extrapolate the 
dust opacity \ v at 1-2 mm from its value at 250 /im, 
i.e. X250^m = 10 cm 2 g _1 , assuming the power-law de- 
pendence Xv °c v v 1 where p is the emissivity index and 
strongly depends on the mineralogical composition of the 
grain and on their physical shape. 

Under the hypothesis of optically thin emission, the 
emissivity index p may be derived from the spectral index 
in the millimetric and submillimetric spectral range, where 
the dust emits as a blackbody modified by the frequency 
dependent dust opacity, that is F v oc XvB u (T d ). 

We thus proceeded iteratively by fitting a modified 
blackbody to the infrared and millimetric data to estimate 
the dust temperature, while the emissivity index has been 
derived from a linear fit in the log v - log B ^f d ) +° ^ ne m " 
tegrated flux density from 100 to 1.2 mm. In the case 
of IRAS 19075+0432, it is to note that the MSX data 
at 8.28 /Ltm suggest the presence of both warm and cold 
circumstellar dust. 

For sources with only an upper limit on IRAS 100 /im 
flux, the same analysis as above was performed using the 



IRAS 60 /im flux to constrain minimum value for p. In 
n g- the resulting curves are overplotted to the observed 
SEDs. 

In table |3| the best fit parameters are listed together 
with the derived dust opacities and masses calculated from 
the eq. mfor A = 1.2 mm. 

3.2. SED modeling with DUSTY code 

For six of the detected post-AGB stars, the observations 
scan a wide spectral range between optical bands and ra- 
dio regions. All the available informations on the sources, 
including the results from SED fits with envelope models, 
have been collected from the literature, in order to verify 
the compatibility of the previously determined parameters 
with our 1.2 mm fluxes. 

The observed SEDs are plotted in fig. |3 along with 
best fit SEDs from the radiative transfer code DUSTY 
(Ivezic et al. I1999fl . The code allows to calculate the 
radiative transfer through a spherically symmetric dust 
shell and determines the spectral energy distribution on 
the basis of specified properties of the radiation source, 
dust composition and dust shell properties. In particular, 
DUSTY allows to use six different types of dust grains 
which are assumed to be distributed in a spherical shell. 
For all the sources we assumed Planckian SED for the 
central star and the modified MRN (Mathis, J.S., Rumpl, 
W., Nordsieck, K.H. 1977) power law for the grain size 
distribution n(a) oc a~ 9 , where a is the grain size and q 
is fixed to 3.5, as it is commonly used, if not otherwise 
specified in the text. For more details about the other op- 
tions we reference to the description in the user manual. 
When available, we calculated the SED by assuming the 
stellar and dust parameters taken from previous references 
in literature, otherwise the best-fit model parameters are 
sought iteratively by fitting the shape of the observed 
SED. With the aim of add more observative constrains, 
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IRAS 07134+1005 
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Fig. 2. Left: Image of IRAS 07134+1005 observed at 11.9 fim by Hony et al. ip00l)l . Right: The map simulated at 
11.9 /Ltm using the DUSTY code, assuming the parameters reported in table The map has been convolved with a 
circular beam having HPBW of 0".83. In both the maps the contours indicate the 5-95% intensity levels in steps of 
10%. 



we compared the infrared maps of the sources available 
in literature with the simulated maps obtained from the 
code at the same observative frequencies. All the sim- 
ulated maps have been convolved with the beam 
of the instrument used to obtain the map used 
for comparison. An example of such a comparison is 
shown in fig-Hin the case of IRAS 07134+1005. In partic- 
ular we refer to Hony et al. l|2UU5|l for IRAS 07134+1005, 
Gledhill et al. for IRAS 17436+5003, Jura 

& Werner i flMty for IRAS 19114+0002, Ueta et 
al. PUUT|> for IRAS 22272+5435 and to 2MASS 
Atlas Image provided from Infrared Science Archive 
(bungo.ipac.caltech.edu/applications/2MASS/IM/ inter- 
active.html) for IRAS 20000+3239 IRAS 20028+3910. We 
thus rejected all those solutions that involve shell sizes 
which don't match the observative constrains according 
to the rms noise level of the measures. 

The remaining four sources of our sample, detected in 
the survey, were not fitted because of the lack of either 
data and knownledge of stellar parameters, involving too 
many degrees of freedom. 

A summary of the adopted input parameters used to 
produce the model fits in fig. 03 is given in table 0| In 
addition, we report the derived stellar and dust shell pa- 
rameters beside the physical quantities used in scaling the 
dusty outputs in accordance with the dusty manual (Ivezic 
et al. I1999fl . In all calculations we assumed a density for 
the grain material of 3 g/cm 3 and a gas/dust ratio of 220. 

The fit's parameter that is mainly constrained from the 
millimetric and submillimctric measures is the shell outer 
radius, as these observations probe the cool, outer parts 
of the CSE. However, the value derived from the SED fit, 
and reported in table 21 bas to be meant as an estimate 
at the distance from the star where the wind density is 
close to the ambient density. Thus, the modeled envelope 



could result to be more extended than the observed one 
if its outer parts comes to have surface brightness lower 
or equal to the rms noise level. In particular, it is not 
surprising that the radii derived for the CSE outer edge 
result to be much greater than those derived from the near 
or mid infrared observations, since at those frequencies 
most of the flux comes from the hotter inner regions of 
the envelope. 

On the basis of the mass loss rate and shell size 
derived from the SED fit, we calculated the total 
envelop mass, assuming a constant wind expansion 
velocity which has been taken from references in 
literature. The dust envelope masses derived for a 
gas/dust ratio of 220 are also reported in table QJ 
Such values are consistent with those reported in 
table|31within a factor of 2. We can thus asses that 
there is a good agreement between the two dust 
mass estimates, expecially given the rough approx- 
imations made in the assumption of constant mass 
loss rate and in the choice of the gas/dust ratio, 
which could differ significatively from the assumed 
value as a function of the C/O ratio. 

The fit of the model to the observed flux distribution 
for each object is discussed in the next section. 

4. Results and discussion 

IRAS 07134+1005. This object is known to have strong 
features near 21 and 30 /iin, whose carrier species have not 
been firmly identified (Kwok et al. H989f) . It is an F5 super- 
giant having C/O « 1 (Van Winckler & Reyniers [2000 ) . 
Despite of the very simplified dust composition used in 
our model calculations, the parameters derived from our 
fit are in good accordance with the ones calculated from 
Meixner et al. 120041 and Hony et al. (2003), using different 
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Fig. 3. Spectral energy distributions of the emission towards six stars in the sample. The asterisks indicate measures 
from this work. The dashed lines show the SEDs computed using the DUSTY code. 



dust radiative transfer programs, with the exception of the 
mass loss rate that is larger for a factor of 3 from the one 
calculated from Hony et al. |2003). On the other hand, 
our derived value for M agrees with the values obtained 
from the SED model performed by Hrivnak et al. (2000) 
and from Jura et al. I|2000|) on the basis of their mid-IR 
images of the source. 

IRAS 17436+5003. It is a F3 lb supergiant with an 
oxygen-rich chemistry (Justtanont et al. ©92). Its SED 
has been modeled by several authors adopting different 
choices mainly in the selection of grain parameters. While 
Hoogzaad et al. I|2002|l and Meixner et al. Q2002fl adopt the 
usual grain size distribution n oc a~ 3 ' 5 with a minimum 
size of 0.18 /im and 0.2 /im respectively, Gledhill & Yates 
( 2003J) opt for very small grains, using a steep size distri- 
bution n oc a~ 6 and minimum grain size 0.01 /im, that, as 
pointed out from the authors, better account for both the 
SED shape within the whole frequency range and the ob- 
served high degrees of near-IR linear polarization. In our 
calculation we followed the latter way. As a consequence, 
the dust mass reported in tab. could be underestimate 
because the assumed dust opacity is derived from an ap- 
proximated grain model with an average grain size of 0.1 
/xm ( Hildebrand, 1983). This value is a factor of 5 greater 
than the grain size obtained from the n oc a~ 6 distribution 



in the adopted range for a. Assuming for such small grain 
radius Xs50^m = 0.54 cm 2 g" 1 (Gledhill & Yates HSH , 
the derived dust mass increases of one magnitude order. 

The inner shell radius derived from the SED fit agrees 
with the value of 1400 AU (2.1 x 10 16 cm at 1.2 kpc) 
obtained by Hoogzaad et al. (|2002[) . but is about 30 - 40% 
larger than the same parameter derived from Gledhill & 
Yates (|2003[1 and Meixner et al. H2002|l . who nevertheless 
used an asymmetric dust model. The outer radius is larger 
than the values obtained from Hoogzaad et al. (|2002|l and 
Meixner et al. <|2002[) . but close to that derived by Gledhill 
& Yates H2003fl and to the extension of 6". 5 derived from 
the CO observation (Bujarrabal et al.,[1992). 
IRAS 19114+0002. It is an oxygen-rich object of spectral 
type G5 la (Hrivnak et al. 1T2S2J), but the real evolutiona- 
ry status of this object is still controversial. It has been 
classified either as post-AGB star having a luminosity of 
about 10 4 L© and lying at 1 kpc (Hrivnak et al. 198§|) or 
as a massive red supergiant with luminosity near 3 x 10 5 
L0 and lying to a distance of 6 kpc (Hawkins et al.[l995). 

The expansion velocity of the circumstcllar envelope 
v c « 35 km s _1 determined from the profile of the circum- 
stcllar CO emission (Zuckerman & Dyck 1986. Bujarrabal 
et al. I1992|) is significantly higher than the typical value 
of 15 km s" 1 for an AGB star and favors the supergiant 
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Table 4. Input and derived stellar and dust parameters resulting from the fit to the SED. 



Property 


IRAS 07134+1005 


17436+5003 


19114+0002 


20000+3239 


20028+3910 


22272+5435 


T*[K] 


7250 


7500 


5660 


5000 


7000 


5300 


Chem. 


C 


O 


O 


C 


C? 


C 


Td[K] 


150 


120 


110 


170 


180 


160 


L [L Q ] 


6000 


3440 


300000 


630 


6600 


6700 


Av 


0.5 


1.2 


2 


2.5 


1.3 


2.5 


r; [cm] 


4.76 x 10 16 


1.91 x 10 16 


1.73 x 10 17 


8.18 x 10 15 


1.88 x 10 1B 


2.46 x 10 16 


To [cm] 


3.48 x 10 17 


1.53 x 10 17 


6.88 x 10 17 


2.05 x 10 17 


1.88 x 10 17 


3.11 x 10 17 


M M^] 

L yr 1 


3.48 x 10" 5 


3.48 x 10~ 5 


1.09 x 10" 3 


6.92 x 10" 6 


8.93 x 10 _s 


3.12 x 10 -5 


v exp [km/s] 


10 .0 


15.5 


35.0 


12.0 


16.0 


10.0 


M d [M Q ] 


1.66 x 10~ 3 


4.35 x 10~ 4 


2.31 x 10~ 2 


1.63 x 10~ 4 


1.35 x 10~ 3 


1.28 x 10~ 3 


Ref. 


1, 2, 3, 4 


5, 6, 7, 18, 19, 20 


3, 8, 9, 10, 18, 19, 21 


11, 12, 13, 19, 22 


14, 15, 18, 19, 20 


13, 16, 17, 18 



References - (1) Meixner et al. 12001 (2) Hony et al. 120031 (3) Hrivnak et al. rTOSOl (4)Van Genderen et al. 1T087)1 (5) Ghedill 
et al. 130031 (6) Ueta et al. 120001 (7) Skinner et al- ITOMl (8) Van der Veen lTTOl (9) Hawkins et al. 1X0051 (10) Thevenin et al. 
120001 (11) Kwok et al.. lT9"9"5l (12) Hrivnak et al. HOOSI (13) Volk et al. 120021 (14) Su et al. I200T1 (15) Bujarrabal et al. I200T1 
(16) Hrivnak et al. lTOOTl (17) Ueta et al. l200Tl (18) Walmsley et al. HOOTl (19) Gledhill et al. 120021 (20) Likkel et al. 116611 
(21) Bujarrabal et al. 119921 (22) Hrivnak et al. 120001 



hypothesis, that is though still doubt (Josselin & Lebre 

As for IRAS 17436+5003, Gledhill & Takami ipMjl 
pointed out the need to adopt a steep power law for the 
grain size distribution, in order to agree with the high de- 
grees of linear polarization observed in the near-IR. We 
have thus adopted the grain size distribution n oc a~ 6 
with a m i n = 0.005 fim and a max = 0.25 /im. The SED 
fit gives a r- m = 1.7 x 10 17 cm, which is in good accor- 
dance with the inner shell radius of about 1.5 x 10 17 
cm obtained from both mid-IR images (Jura & 
Werner 1999) and near-IR imaging polarimetry 
observations (Gledhill et al. 2001). Such value is 
very close to the inner shell size of 1.8 x 10 17 cm 
derived, for the assumed distance of 6 kpc, from 
high resolution CO observations performed from 
Jura et al. (200lJ. 

The dust mass obtained is in good accordance with 
the one derived from Gledhill et al. (2002) on the basis of 
submillimeter observations. 

IRAS 20000+3239. Low resolution K-band spectra per- 
formed by Davis et al. I}2003}1 showed a rather compact 
source with angular size lower than 2", as well Hrivnak 
et al. H1999J1 measured 1".6 in diameter in V band. 
Comparing with the SED parameters derived from Volk et 
al. H2002|l . we obtain good accordance for both the mass 
loss rate and inner shell radius, as well as for the shell 
mass derived from Gledhill et al. |2002). 
IRAS 20028+3910. This object is characterized from a 
bipolar morphology, with the central object highly ob- 
scured in optical and near- IR (Su et al. I3TM1 Ueta et 
al. I2000[> . The SED, constructed with the data reported 
by Su et al. (20QIJ and references therein, thus shows a 
mid and far infrared peak much most bright than the one 
in the near infrared. Neri et al. H1998J1 fitted the CO 1-0 
visibility data with an elliptical gaussian component with 



a size of 3.5" x 11.1". The dust mass obtained from our 1.2 
mm measures agrees very well with the value calculated 
from previous submillimeter data (Gledhill et al. 2002). 
IRAS 22272+5435. It is an extremely carbon rich object 
and, as IRAS 07134+1005, shows peculiar infrared spec- 
tral features. From the analysis of CO 1-0 visibility data, 
Neri et al. (1998) measured an extended envelope size 
(FWHM) of 21". From subarcsecond mid infrared ima- 
ging study (Ueta et al. 12001}! the dust shell was found to 
have a toroidal structure with a 0.5" inner radius, which 
corresponds to l.lxlO 16 cm at 1.6 kpc. Our fit indicates 
an inner radius which is consistent with that measured by 
Ueta et al. (.2001 ) to whithin a factor of 2. A major dis- 
crepancy is found in the estimates of mass loss rate. On 
the basis of their radiative transfer calculations, in fact, 
Ueta et al. (|2001}l derived a wind with M = 4.1x 10~ 6 
M© yr -1 . Our value is closer with the mean mass loss 
rate obtained from the CO study performed Bujarrabal et 
al. (2001), which is 1.8 x 10~ 5 M© yr -1 , when scaled for 
our assumed distance. The authors report a total nebular 
mass of 0.18 M©, that gives a dust mass of 8.0 x 10~ 4 M© 
for an assumed gas to dust ratio of 220, and which is in 
good agreement with our estimate. 

5. Summary and outlook 

We have presented the results of 1.2 mm continuum obser- 
vations for a sample of 24 sources classified as post-AGB. 
Continuum emission was detected toward 11 objects, while 
uncertain detection is reported for two other sources. 

The circumstellar dust masses were derived from our 
1.2 mm measures, assuming that the emission is due to 
optical thin dust. For the sources of which the distance is 
known, the circustellar dust masses range between about 
6 x 10~ 4 M© and 2.4 x 10 -3 M©, with exception of 
IRAS 19114+0002, whose post-AGB nature is, however, 
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still in question. For the other objects we derived lower 
dust masses, indicating a more young circumstellar enve- 
lope, but the errors in the choice of the common assumed 
distance could have effected the derived values. 

For six of the detected sources, we compared the ob- 
served SEDs, constructed with additional data from lite- 
rature, and the model spectra obtained using a symme- 
tric radiative transfer code. This allowed us to estimate 
some physical parameters of stars and envelopes which 
have been compared with previously results reported in 
literature. 

The high detection rate (» 46%) seems to support the 
goodness of our selection rules and an extension of the 
millimetric survey to the remain targets in our sample 
is necessary. Stars in our full sample belong to different 
evolutionary phases in the transition from AGB to PN, as 
the stars with optical counterpart should be more evolved. 
Once the full sample will be observed, the comparison be- 
tween the derived physical properties of different envelopes 
will provide fundamental informations on this evolutiona- 
ry phase not fully understood yet. 

We still note that the CSEs surrounding such kind of 
objects appear as good targets for the first light projects 
for future millimeter arrays such as the Atacama Large 
Millimeter Array (ALMA). From our fits we derive typical 
dimensions of CSEs ranging from 1.5 x 10 17 to 1.2 x 10 18 
cm, which, combined with the distances as reported in ta- 
ble 3, correspond to angular sizes from few arcsecs up to 
« 13". This implies that such CSEs can be, in principle, 
mapped, in great details with the foreseen ALMA angu- 
lar resolutions by using a combination of both compact 
(to fully recover all the flux) and extended configurations 
(Wilson et al. 2005). The foreseen capabilities of ALMA 
will allow, at least for the more compact sources of our 
sample, and in general for post-AGB sources, to directly 
map the dusty envelopes at several millimetric and submil- 
limetric frequencies. This could provide a better constrains 
to the modeling of CSE. Furthermore, a detailed map of 
CSE could evidence any kind of structured morphologies 
that can be related to different mass loss episodes, suffered 
by the star during the AGB evolutionary phase. 

To evaluate the possibility to actually resolve and 
map the CSEs of our sample, we need to compare, 
at each frequency channel, the expected flux densities 
with the foreseen ALMA sensitivity. From the fitted 
SEDs we have then extrapolated the flux densities in 
the ALMA first light channels, namely 0.5, 0.6 0.9 
1.3, 2 and 3 mm. Expected ALMA sensitivities have 
been calculated by using the ALMA Sensitivity Calculator 
(www.eso.org/projects/alma/science/bin/sensitivity.html) 
in the case of first- light, assuming that only 8 antennae 
will be available, and in the case of the full (64 antennae) 
array. For both configurations we derived a detection rate 
close to 100% over the 3cr at almost all frequencies. In 
fig. 0] we show the expected percentage of the studied 
objects that ALMA will allow to observe with dynamical 
range greater than 50. It is evident how ALMA will allow 
us not only to better sample the millimetric range of the 
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Fig. 4. Percentage of objects, detected at 1.2 mm, which 
are observable with dynamical range grater than 50, as- 
suming calculate ALMA sensitivity with both 8 and 64 
antennae. 

source SEDs, that up to now is very poor, but also to 
obtain, in most cases, multifrequencies high resolution 
maps of the circumstellar matter surrounding the stars 
of our sample, extended to the sources belonging to the 
south hemisphere. 
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prove this paper. 

References 

Balick, B., Frank, A. 2002, ARA&A, 40,439 

Bujarrabal, V., Alcolea, J., Planesas, P. 1992, A&A, 257, 701 

Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., Sanchez 

Contreras, C. 2001, A&A, 377, 868 
Clark, J. S., Egan, M. P., Crowther, P. A, Mizuno, D. R., 

Larionov, V. M.,Arkharov, A. 2003, A&A, 412, 185 
Davis, C. J., Smith, M. D., Stern, L., Kerr, T. H., Chiar, J. E. 

2003, MNRAS, 344, 262 
Frank, A., Balick, B., Livio, M. 1997, Planetary nebulae, IAU 

Symp. 180, Habing, H.J. and Lamers H. J. G. L. M. eds., 

225 

Garcia-Lario, P., Manchado, A., Pych, W., Pottash, S.R., 1997, 

A&AS, 126, 479 
Garcia- Segura, G., Lopez, J. A., Franco, J. 2005, ApJ, 618, 919 
Gledhill, T. M., Takami, M. 2001, MNRAS, 328, 266 
Gledhill, T. M., Bains, I., Yates, J. A. 2002, MNRAS, 332, L55 
Gledhill, T. M., Yates, J. A. 2003, MNRAS, 343,880 
Hawkins, G. W., Skinner, C. J., Meixner, M. M., Jernigan, J. 

G., Arens, J. F., Keto, E., Graham, J. R. 1995, 452, 314 
Hildebrand, R. H. 1983, QJRAS, 24, 267 
Hrivnak, B. J., Kwok, S., Volk, K. m. 1989, ApJ, 346, 265 
Hrivnak, B. J., Kwok, S. 1991, ApJ, 371, 631 
Hrivnak, B. J. 1995, ApJ, 438, 341 

Hrivnak, B. J., Langill, P. P., Su, K. Y. L., Kwok, S. 1999, 
ApJ, 513, 421 

Hrivnak, B. J., Volk, K., Kwok, S. 2000, ApJ, 535, 275 



10 



Buemi et al: A 1.2 mm MAMBO survey of Post-AGB stars 



Hrivnak, B. J., Bieging, J. H. 2005, ApJ, 624, 331 

Hony, S., Tielens, A. G. G. M,; Waters, L. B. F. M., de Koter, 

A. 2003, A&A, 402, 211 
Hoogzaad, S. N., Molster, F. J., Dominik, C, Waters, L. B. F. 

M., Barlow, M. J., de Koter, A. 2002, A&A, 389, 547 
Huggins, P. J., Muthu, C, Bachiller, R., Forveille, T., Cox, P. 

2004, A&A, 414, 581. 
Ivezic, Z., Nenkova, M., Elitzur, M. 1999, User Manual for 

DUSTY, University of Kentucky Internal Report 
Josselin, E.; Lebre, A. 2001, A&A, 367, 826 
Jura, M., Werner, M. W. 1999, ApJ, 525, L113 
Jura, M., Chen, C, Werner, M. 2000, ApJ, 544, L141 
Jura, M., Velusamy, T., Werner, M., W. 2001, ApJ, 556, 408 
Justtanont, K., Barlow, M. J., Skinner, C, J., Tielens, A. G. 

G. M. 1992, ApJ, 392,L75 
Kelly, D. M., Hrivnak, B. J. 2005, ApJ, 629, 1040 
Kreysa, E., Gemuend, H.-P., Gromke, J., et al. 1998, SPIE, 

3357, 319 

Kwok, S.,Volk, K. M., Hrivnak, B.J. 1989, ApJ, 345, L51 
Kwok, S., Hrivnak, B. J., Geballe, T. R. 1995, ApJ, 454, 394 
Likkel, L., Forveille, T., Omont, A., Morris, M., 1991, A&A, 
246, 153 

Mathis, J. S., Rumpl, W., Nordsieck, K. H. 1977, ApJ, 217, 
425 

Meixner, M., Ueta, T., Bobrowsky, M., Speck, A. 2002, ApJ, 
571, 936 

Meixner, M., Zalucha, A., Ueta, T., Fong, D., Justtanont, K. 

2004, ApJ, 614, 371 
Neri, R., Kahane, C, Lucas, R., Bujarrabal, V., Loup, C. 1998, 

A&AS, 130, 1 

Oudmaijer, R.D., van der Veen, W.E.C.J., Waters, L.B.F.M., 
Trams, N.R., Waelkens, C, Engelsman, E. 1992 A&A 96, 
625 

Oudmaijer, R.D., 1996, A&A 306, 823 

Pasquali, A., Comeron, F. 2002, A&A, 389, 874 

Pollack, J.B., Hollenbach, D., Beckwith, S., Simmonelli, D. P., 

Roush, T., Fong, W. 1994, ApJ, 421, 615 
Pottasch, S.R., Parthasarathy, M. 1988, A&A, 192, 182 
Sahai, R., Trauger, J. T. 1998, AJ, 116, 1357 
Sahai, R., 20000, Post-AGB Objects as a Phase of Stellar 

Evolution, R. Szczerba and S. K. Grny eds., Ap.S.S.L., 265, 

53 

Skinner, C. J., Meixner, M. M., Hawkins, G. W., Keto, E. 

,Jernigan, J. G., Arens, J. F. 1994, ApJ, 423, L135 
Su, K. Y. L., Hrivnak, B. J., Kwok, S. 2001, AJ, 122, 1525 
Thevenin, F., Parthasarathy, M., Jasniewicz, G. 2000, A&A, 

359, 138 

Trams N.R., Waters, L.B.F.M., Waelkens C. et al., 1991, A&A 
SS 87, 361 

Ueta, T., Meixner, M., Bobrowsky, M. 2000, ApJ, 528, 861 
Ueta, T., Meixner, M., Hinz, P. M., Hoffmann, W. F., 

Brandner, W., Dayal, A., Deutsch, L. K., Fazio, G. G., 

Hora, J. L., 2001, ApJ, 557, 831 
Umana, G., Buemi, C. S., Trigilio, C, Leto, P. 2005,A&A, 437, 

LI 

van der Veen, W.E.C.J., Waters, L.B.F.M., Trams, N.R., 

Matthews, H.E., 1994, A&A 285, 551 
van Genderen, A. M., Greidanus, H., van Driel, W. 1986, A&A, 

155, 72 

van Winckel 1997, A&A, 319, 561 
van Winckel, H., Reyniers, M. 2000, A&A, 354, 135 
Volk, K., Kwok, S., Hrivnak, B. J., Szczerba, R. 2002, ApJ, 
567, 412 



ilson, T. L., Beasley, A. J., Wootten, H. A. , 2005 in The Cool 
Universe: Observing Cosmic Dawn, ASP Conference Series, 
Vol. 344, 232 

Walmsley, CM., Chini, R., Kreysa, E., Steppe, H., Forveille, 

T., Omont, A., 1991 A&A, 248, 555 
Zuckerman, B., Dyck, H. M. 1986, ApJ, 311, 345 

Zyl ka, R. 1998, Pocket Cookbook for MOSPI Software, 

http:/ /iram.fr/IRAMES/otherDocuments/manuals/Datared/pockcoo.ps 



